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ABSTRACT 


A  formalism  is  presented  which  describes  the  response  of  a  complex  of  coupled  one¬ 
dimensional  dynamic  systems  to  an  impulse  drive.  The  formalism  is  based  on  an  impulse  response 
operator  which  relates  a  drive  applied  to  one  point  in  the  complex  to  the  response  at  any  point  in  the 
complex.  The  formalism  is  derived  directly  in  the  time  domain  and  the  impulsive  drives  which  can 
be  accommodated  must  be  finite  in  time  and  applied  at  a  spatial  point.  The  constituent  systems 
must  be  one-dimensional  and  possess  a  pulse  propagation  velocity  which  is  not  a  function  of 
position  within  the  system.  Systems  interact  through  junctions  which  also  define  their  spatial 
extents.  The  junctions  are  characterized  by  reflection  and  transmission  coefficients  which  modulate 
the  amplitude  of  reverberant  components  and  by  delays  in  the  reflections  and  transmissions. 
Propagation  in  the  systems  is  characterized  by  losses.  Several  simplistic  examples  are  calculated 
and  presented  to  illustrate  the  type  of  information  which  the  formalism  can  provide. 


ADMINISTRATIVE  INFORMATION 

This  work  was  supported  by  the  Propulsion  and  Auxiliary  Systems  Department,  Code  27, 
of  the  David  Taylor  Research  Center. 


INTRODUCTION 

This  paper  introduces  a  formalism  which  describes  the  temporal  response  of  ?  complex 
consisting  of  multiple  connected,  one-dimensional  systems.  The  formalism  is  based  on  a  matrix 
impulse  response  operator 

g(x|x'.t|  O  =  (gji(Xj  |Xj',  t  It'))  ^ 

which  depends  solely  on  the  properties  of  the  complex.  In  particular,  g  depends  on  the 
propagation  properties  of  each  of  the  systems,  the  boundary  conditions  for  each  system,  the  spatial 
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coordinate  and  time  of  the  observation  and  drive  in  each  system.  It  is  assumed  that  the  impedance 
operator  governing  the  system  is  linear  so  that  the  response  of  the  complex  p  (x,  t),at  some  point  x 
and  at  time  t,  can  be  given  by  the  operation  of  g  on  an  external  drive,  (i.e.  one  which  does  not 
interfere  with  or  depend  on  the  response  of  the  complex),  in  the  usual  way: 

P(x,  t)  =  J  gOc  |  ,  t  |  t')  d&'  dt'  geU' ,  t')  , 

<&'  =  dx,r  «ij  •  (2) 

The  rank  of  the  matrix  and  vector  quantities  in  equation  (2)  equals  the  number  of  systems  in  the 
complex,  and  individual  components  of  p  ,  say  Pj  (Xj,t),  describe  the  response  in  the  (i)th  system  at 
the  point  Xj  and  at  time  t.  Similarly,  each  component  of  the  drive,  e.g.  pej  (xj,  t)  describes  the 
drive  in  one  system,  e.g.  the  (j)th,  at  a  point,  e.g.  Xj,  and  at  time  t.  Each  component  of  the  vector 
spatial  variables  x  and  x'  is  the  spatial  positions  of  observation  and  drive  respectively  in  each  of 
the  systems.  The  response  due  to  a  drive  in  the  (j)th  system  at  the  point  xj  and  at  time  t'  is  given  by 
the  element  gjj  (xj  |xj,  t  |t'),  the  transfer  impulse  response  operator,  which  makes  the  connection 
and  completely  accounts  for  the  dynamic  interaction  of  all  members  of  the  complex. 

The  focus  here  is  on  the  temporal  response  of  the  complex,  and  the  development  of  the 
appropriate  impulse  response  operator  will  closely  parallel  the  derivation  of  the  spatial  impulse 
response  function  [1,2].  This  approach  tracks  a  single  spectral  component  in  the  response  as  this 
component  interacts  with  the  complex  and  superposes  the  results  of  these  interactions  at  an 
observation  point  One  of  the  useful  characteristics  which  the  present  formalism  shares  with  this 
approach  is  that  the  impulse  response  operator,  and  hence  the  response,  so  derived,  is  obtained  as 
a  sum  of  two  direct  terms  and  four  reverberant  terms.  The  reverberant  terms  express  the  response 
arriving  at  the  observation  point  from  each  direction,  and  in  general,  each  of  these  are  the 
superposition  of  two  components  which  left  the  source  in  each  direction.  This  division  into 
components  can  be  shown  by  writing  equation  (2)  in  the  expanded  form: 
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(3a) 


PaU,t)  =  Jg«OUi',t|t')  d&'dt'k(&',  t')  , 

p(x,  t)  =  {pj(Xj,  t)},  Pe(x',t')  =  {peiCx'i,  t')}  , 

g&lx' .  1 1 1')  =  (gji(xj  |  x; ,  1 1 1'))  , 

where  a  =  r  or  q,  r  and  q  indicating  opposite  ends  or  directions  in  a  system.  Thus  for  example, 

Pj  (xj  *  0  is  the  response  in  the  (j)th  system  which  is  traveling  towards  the  r  junction.  The  external 
drive  in  each  system  is  composed  of  two  components,  one  launched  in  each  direction, 


$e  —  ^  CJ  Pea  • 
~  a=r,q  ~  ~ 


(3b) 


The  two  components  of  the  impulse  response  operator  in  equation  (3a)  are  further  composed  of 
components  which  relate  drive  components  emanating  in  each  direction  to  the  arrival  components; 
specifically. 


(3c) 


where,  if  a  =  r  then  P  =  q  and  vice-versa,  and  g“  is  the  direct  component  of  response;  i.e.  the 
response  which  has  not  interacted  with  the  boundaries.  Also,  for  example,  the  gg  component  of 
the  impulse  response  function  relates  the  p  component  of  the  drive  with  the  a  component  of  the 
response. 


When  formulating  the  response  in  this  way,  it  is  often  assumed  that  the  complex  is 
stationary  in  time.  Under  this  condition  one  can  write 

g(x|x',t|t')->  (2u)-1/2  g(x|x',  t - 1')  . 


(4) 


Equation  (2)  can  then  be  Fourier  transformed  in  the  temporal  variable  yielding 

p(x,  co)  =\ g(x  t  x',co)dx,pe(x/,  co) 


(5) 
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where  the  frequency,  co,  is  the  Fourier  conjugate  of  the  temporal  variable,  t,  and 
p(x,co),pe(x,  co)  and  g(x  j  x' ,  co)  5(co-  co')  are  the  Fourier  transforms  of  p(x ,  t),  pe(x,  t) 
and  ( 2k )'1/2  g(x  |  x' ,  t  - 1')  respectively.  It  is  then  observed  that  under  these  conditions  of 
temporal  stationary,  equation  (2)  can  be  transformed  to  an  equation  which  is  algebraic  in  the 
frequency.  Without  this  stationarity,  equation  (5)  must  remain  operational  in  the  frequency 
variable.  In  the  present  development,  the  impulse  response  operator  will  be  manipulated  in 
operator  form  and  the  complex  will  not  be  considered  to  be  stationary  in  time.  The  Fourier 
transformation  in  that  variable  will  not  be  performed;  rather  we  will  consider  the  nature  of  g  in  the 
temporal  domain  with  the  drive  points  and  the  observation  points  held  constant.  The  impedance 
operator  for  the  complex  is  linear;  thus  the  response  to  a  spatially  distributed  source  can  be 
constructed  as  a  summation  (integration)  in  the  spatial  domain  in  an  analogous  way  in  which  the 
Fourier  components  are  conventionally  summed  (integrated)  to  get  the  response  to  a  drive 
distributed  in  time. 

Equation  (2)  will  hereafter  be  written  as 

£(&  0  =  fe(x,  x'»  t  |t')  pe(x',  t')  ,  (6a) 

with  h  being  termed  the  impulse  response  operator.  The  major  object  of  this  paper  is  to  present  an 
explicit  form  for  h  and  to  explore  some  of  its  characteristics  in  the  temporal  domain  by  picking 
Fixed  drive  and  observation  points,  thus  effectively  "putting  a  hold"  on  its  spatial  characteristics.  A 
comprehensive  formalism  which  treats  the  spatial  and  temporal  aspects  in  a  unified  manner  is  in 
preparation  [3]. 

The  temporal  impulse  response  operator  for  a  single  system  would  be  the  same  as  equation 
(5)  but  the  matrix  and  vector  quantities  would  all  be  of  rank  one,  i.e.  for  a  single  system, 

p(x,  t)  =  h(x,  x',  t  |t')  pe(x',  t')  .  (6b) 
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Since  the  response  to  a  temporal  pulse  is  more  easily  visualized  in  a  single  system  than  in  a 
complex,  the  formalism  will  first  be  developed  and  demonstrated  by  calculation  for  a  single 
system.  After  this,  the  extension  to  multiple  systems  is  straightforward;  this  will  be  done  and 
demonstrated  in  the  final  two  sections.  There  are  several  assumptions  about  the  systems  which  are 
basic  to  both  the  single  and  multiple  system  formulation.  Specifically,  the  impedance  operator 
which  governs  the  system  is  assumed  to  be  linear  so  that  solutions  may  be  superposed,  and  it  is 
assumed  that  the  systems  are  one-dimensional;  i.e.  if  the  physical  system  being  modeled  is 
spatially  two  or  three  dimensional,  then  the  impedance  operator  for  the  system  must  be  separable  in 
such  a  way  that  the  spatial  dimension  of  interest  can  be  treated  independent  of  the  others.  Also,  a 
system  can  have  only  one  propagation  speed  and  this  must  be  constant  throughout  the  spatial  extent 
of  the  system.  Separate  systems  can,  of  course,  possess  different  speeds,  and  within  a  particular 
system,  the  speed  in  one  direction  may  differ  from  the  speed  in  the  other  direction.  In  view  of  this 
last  assumption,  a  physical  system  with  more  than  one  mode  of  propagation,  and  therefore  more 
than  one  type  of  pulse  propagation,  would  be  formulated  as  several  systems,  one  for  each  mode, 
which  are  allowed  to  interact  with  each  other  only  at  the  junctions. 


THE  DELAY  OPERATOR 

The  temporal  impulse  response  operator,  h  in  equation  (6),  uses  the  delay  associated  with 
the  propagation  of  a  pulse  between  two  spatial  points,  x'  and  x,  as  the  "building  block"  in  an 
analogous  way  in  which  the  spatial  formalism  [1,2]  uses  spatial  propagators.  The  spatial 
formalism  defined  a  spatial  propagator  in  terms  of  the  eigenfunctions,  <Kx),  for  a  specific  response 
mode  of  a  spatially  one-dimensional  system  with  no  boundaries.  This  propagator  is  of  the 
form;  [4] 

t(x|x')  =  Wx')'1  <J>(x)  •  (7) 
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For  purposes  of  constructing  the  temporal  impulse  response  operator,  we  consider  the  delay 
experienced  by  the  propagation  of  a  pulse  in  a  system  which  has  the  same  propagation 
characteristics  as  the  system  of  interest  but  which  has  no  boundaries;  i.e.  a  specific  response  mode 
in  a  spatially  one-dimensional  system  which  is  infinite  in  extent.  If  the  propagation  speed  of  this 
pulse  is  u,  then  the  delay  in  propagating  from  the  point  x'  to  x  is 

T(x|x')  =  |x'-x|/u  .  (g) 

We  define  a  "delay  operator",  A(t),  such  that 

A(t)  f(t)  =  f(t  +  x), 

A(ti)  A(x2)  =  A(ti  +t2)  •  (9) 

The  following  sections  will  use  the  above  delay  operator  to  construct  the  temporal  impulse 
response  operator,  first  for  a  single  system,  and  then  for  a  complex  of  connected  systems. 


THE  SINGLE  SYSTEM 

A  single  dynamic  system  is  depicted  in  Figure  1.  The  system  is  defined  to  have  a  single 
mode  of  propagation  and  a  velocity  u  associated  with  the  propagation  of  a  pulse;  this  velocity  is  not 
a  function  of  the  spatial  coordinate  x.  In  general,  there  may  exist  different  velocities  associated 
with  the  two  different  directions  of  propagation  within  the  system,  i.e.  ur  and  u<t  with  ur  *  where 
the  superscript  on  u  indicates  the  direction  towards  which  the  pulse  propagates,  [cf.  Figure  1 .] 

The  velocity  u«,  a  =  r  or  q,  is  complex,  ua  =  u$  (1  -  iqa)  where  ria  plays  the  role  of  a  "loss 
factor"  associated  with  propagation.1 

'For  a  system  which  is  non-dispersive,  k  =  co/cp  where  cp  is  a  phase  velocity  independent  of  co,  then 
u  =  cp,  and  setting  u  =  Uq  (1  -  iq)'1  is  equivalent  to  the  more  conventional  practice  of  incorporating  a 
loss  factor  in  the  wavenumber,  k  =  Iq,  (1  -  iq).  When  dispersion  is  present,  the  equivalence  breaks 
down  and  the  loss  factor  associated  with  the  group  velocity  would  not  be  the  same  as  the  one 
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The  end  points  of  the  system  are  determined  by  the  coordinate  values  {xq,  xr},  hence  the 
length  of  the  system  is  L  =  |Xj.  -  xq|  and  the  time  for  propagating  a  pulse  over  the  entire  system  is 
x  (xr|Xq)  =  L/ur  in  the  q  — >  r  direction  and  x  (xq|xr)  =  L/uQ  in  the  r  — >  q  direction.  The  end  points 
or  boundaries  of  the  system  are  characterized  by  reflection  coefficients,  A,,  and  and  by  delays 
experienced  by  an  incident  pulse  daring  reflection,  \  and  xq.  Both  the  reflection  coefficients  and 
the  junction  (reflection)  delays  may  be  complex. 

An  impulsive  drive,  i.e.  one  finite  in  time,  is  applied  to  the  system  at  x'.  In  the  temporal 
coordinate,  the  impulse  is  centered  about  the  point  t'.  The  response  of  the  system  may  be  viewed 
as  a  spatially  finite  wave  packet  whose  "center  of  gravity"  propagates  away  from  x',  in  both 
directions,  with  velocities  ur  and  u<i.  The  point  of  view  adopted  for  the  current  development  is  to 
pick  an  observation  point  x  and  derive  the  response  as  a  function  of  time;  i.e.,  we  seek  to  derive 
the  operator  h  of  equation  (6). 

The  impedance  operator  for  the  system  is  assumed  to  be  linear  and  the  response  is  sought 
as  a  superposition  of  a  "direct"  response  and  a  "reverberant"  response,  with  a  corresponding 
division  of  the  operator  h, 

h(x,x',t|t')  =  h<i(x|x',  t  |t')  +  hrev(x|x',  1 1 1')  .  (10) 

The  "direct"  response  is  given  by  using  hd  in  equation  (7)  and  is  the  response  which  emanates  from 
the  drive  at  {x',  t'}  and  propagates,  without  boundary  interaction,  to  {x,t}.  Using  the  delay 
operator  (1 1),  the  direct  response  is, 

p(x,t)  =  hd(x,x',  t|t')  pe(x', t') 

=  A[x(x  |x')]  Per(x',  t')  U(x-x')  +  A[x(x|x')]  Peq(x',  t')  U(x'-x)  ,  (Ua) 


associated  with  the  wavenumber  and  would  need  to  be  determined  based  on  the  dispersion 
characteristics  of  the  system.  For  purposes  of  illustrating  the  formalism,  a  simple  loss  factor,  Tj,  is 
associated  with  u. 
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where  the  subscripts  r  and  q  on  the  drive,  p^x',  0.  differentiate  between  drive  components  which 
emanate  from  the  point  x'  towards  the  r  and  q  junctions  respectively.  In  general,  pr  *  pq.  Also  in 
equation  (11a),  it  is  assumed  that  xr  >  Xq  and  the  step  function, 


U(x'-x) 


I  1,  x'  >  x 
\  0,  x'<x 


(Hb) 


is  used  to  select  the  drive  component  which  is  operative  in  the  direct  response. 

The  reverberant  contributions  in  equation  (10)  can  be  built  up  by  referring  to  Figure  (1)  and 
noting  that  the  response  at  the  point  x  which  has  been  generated  by  the  r-  component  of  the  drive 
and  which  has  interacted  once  with  the  r-  junction  is 

Pr  (x,  t)  =  A  [x  (x  |  Xr)]  Af  A(xr)  A[x(xr|x')]  Per(x',  t')  .  (12) 


The  superscript,  q,  on  the  response  indicates  that  the  response  is  propagating  toward  the  q  junction 
and  the  subscript,  r,  indicates  that  the  response  was  generated  by  the  r-  component  of  the  drive.  In 
general,  there  is  an  infinite  family  of  such  responses  representing  successive  reverberations  whose 
arrival  times  have  been  delayed  by  successive  complete  "reverberation  times"  for  the  system  and 
have  been  attenuated  by  a  propagation  loss  and  successive  reflections  at  the  boundaries.  The 
reverberation  time  for  the  system  is, 

Xsys  =  Xq  +  Xqr  +  Xr  +  Xrq  , 

Xrq  =  X  (xr|Xq)  =  L/ur  ,  Xqr  =  X  (Xq|xr)  =  L/u^  .  (13a) 


A  reverberation  operator  is  defmed,  which  accounts  for  the  successive  reverberation  delays  and  the 
total  propagation  loss,  and  the  successive  losses  at  the  boundaries: 


2)  —  1  +Ar  Aq  A  (X$ys)  +  [Ar  Aq  A  (Xjys)]  +  [...]  +  ■  •  • 

=  £  [Ar  Aq  A  (XSyS)j 
n  =  0 


(13b) 
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Thus,  equation  (12)  [which  represents  only  the  first  return  of  the  drive  component  which  leaves  the 
point  x'  towards  the  r-junction,  reflects  at  this  junction,  and  goes  to  the  point  x]  can  be  modified  to 
incorporate  all  the  reverberations;  i.e. 

Pr  (x,  t)  =  A[x(x|xr)]  £>  Ar  A(Xr)  A[x(xr|x')]  per(x',  t')  .  (14a) 

Equation  (14a)  is  one  of  the  four  components  in  the  reverberant  response.  In  general,  one  can 
formulate  a  distinct  response  for  each  arrival  direction  and  each  of  these  can  again  be  divided  into 
components  which  have  emanated  from  the  drive  in  the  two  directions.  The  remaining  three 
components  of  the  reverberant  response  can  be  obtained  in  a  similar  way.  They  are: 

Pq(x>  t)  =  A[x(x|xr)]  Ar  A(Xr)  A[x(xr|Xq)]  2)  Aq  A(tq)  A[x(Xq|x)]  Peq  (x  ,  t  )  q4Jj) 


Pr  (x ,  t)  =  A  [x  (x  |  Xq)]  Aq  A  (Xq)  A  [x(Xq  |  Xr)j  2>  Ar  A  (Xr)  A[x(xr|x')]  ^(x'.t')  (14c) 


Pq  (x ,  t )  =  A  [x  (x  |  Xq)]  £)  Aq  A(Xq)  A[x(Xq|x')]  Peq(x',  t')  . 

The  formalism  for  a  single  system  is  now  summarized: 

p(x,  t)  =  h(x,  x',  t ! t')  pe(x%  t')  , 
pe(x',  t')  =  Per(x',  t')  +  peq(x\  t')  ; 

h  (x,  x',  t  1 1')  =  hd  (x,  x',  t  1 1')  +  hrev  (x,  x',  1 1 1')  , 
hd  (x,  x',  1 1 1')  =  A  [x  (x|x')]  =  A  (|x  -  x'|/ur)  U  (x'  -  x) 

+  A  (Ix-x'l/u'i)  U  (x' -  x)  , 

(for  Xr  >  Xq) 
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hrev  -  hf  4-  hq  +  hj  +  h3  =  Z  , 

a,p 

hp  (x,  x',  t  |t')  =  A  [x  (x  I  xa)]  Aa  A(xa)  A  [x  (xa  |  xp)] 

£>ApA(Xp)  A  [x  (xp|x')] 

h«(x,  x',  t  |t')  =  A  [x  (x|xa)]  £)AaA(Xa)A  A[x(xa|x')]  ; 

&  =  n  ^  0  (Ar  A  (^sys)]  > 

^sys  =  +  L/ua  +  Xp  4-  L/liP  , 

L  =  |xr  —  Xq|  . 


RESPONSE  OF  SINGLE  SYSTEM:  AN  EXAMPLE 

The  formalism  presented  in  the  previous  section  is  now  illustrated  by  defining  the 
parameters  for  a  representative  (and  simple)  system  and  performing  two  sets  of  calculations  of  the 
temporal  response  for  two  examples  of  transient  excitations.2  One  excitation  is  a  pulse  which  is 
short  compared  with  the  propagation  time  associated  with  the  system  length,  x0,  and  the  other  is  a 
tone  burst  which  is  relatively  long  compared  with  x0. 

The  system  parameters  which  are  used  for  both  examples  are: 


2The  calculations  were  performed  on  a  desk- top  computer  using  the  symbolic  programming  language 
"Mathematica" .  A  symbolic  processor  of  this  nature  offers  considerable  advantage  in  using  the 
operator  formalism  defined  here. 
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ur  =  =  u, 

Xq  =  0, 

L/u  =  To, 

Tlr  =  TJq  =  TJ  =  0.01, 
Ar  =  0.8,  Aq  =  0.9, 
tr  =  Xq  =  0.0, 


(the  system  is  isotropic), 

(the  temporal  "length”  of  the  system), 
(loss  factors), 

(reflection  coeffecients), 

(no  delays  during  reflection) 


x'/L  =  0.33, 


x/L  =  0.71. 


(16) 


The  drive  used  for  the  first  example  is3 

pe(x',  t')  =  exp  (-20  1 1'| )  cos (50?)  8(x-x')  , 

t'  =  tVt0  ,  (17) 


and  is  shown  plotted  as  a  function  of  time  in  Figure  (2a).  The  time  scale  in  this  and  all  subsequent 
figures  in  this  section  is  normalized  by  %  The  direct  temporal  response  at  the  point  x  (c.f. 
equation  (1 1).)  is  shown  in  Figure  (2b).  The  time  scale  is  condensed  compared  with  that  of  Figure 
(2a)  and  the  origin  (t  =  0)  represents  the  time  t  =  t'  =  0  when  the  center  of  the  pulse  is  applied.  It 
can  be  observed  that  the  pulse  at  x  is  attenuated  compared  with  the  drive  at  x'.  One  of  the 
reverberant  response  components,  p£  (c.f.  equation  (14).)  is  shown  in  Figure  (2c)  on  a  time  scale 
which  is  still  further  condensed  to  show  several  reverberations.  The  decrease  in  amplitude  with 
reflection  and  propagation  loss  and  the  strict  periodicity  in  the  reverberant  response  are  clearly 
evident.  Finally,  the  combined  response  of  the  direct  and  all  four  reverberant  components  is 
shown  in  Figure  (2d).  During  the  first  few  reverberation  times,  the  various  individual  returns  are 
clearly  discemable  and  one  could  identify  an  individual  peak  as  being  a  member  of  a  particular 
family  defined  by  equation  (14).  After  several  reverberation  times  have  elapsed,  the  system  settles 
into  a  modal  response.  This  is  more  clearly  illustrated  in  Figure  (2e)  where  the  total  response  is 


3 The  formalism  does  not  require  that  the  drive  be  a  "good"  function  of  time;  only  that  it  be  represented 
functionally  in  such  a  way  that  the  impulse  response  operator,  h,  sees  the  variable  t  on  which  to 
operate,  and  the  result  be  calculable.  The  drive  could,  for  example,  be  afunctional  representation  of 
an  experimentally  measured  drive. 
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again  shown  on  a  time  scale  which  has  been  still  further  extended  and  the  amplitude  of  the 
response  has  been  plotted  on  a  logarithmic  scale  to  better  accommodate  the  extreme  attenuation. 

For  the  second  example,  a  relatively  long  "tone  burst"  is  used  as  an  external  drive,  and  all 
the  system  parameters  have  been  kept  the  same,  including  the  location  of  the  drive  and  observation 
points.  The  drive  used  is 

pe(x',  t')  =  {exp(— !  t'|)  cos  (20 1')/ (|t' |  +  1)}  S(x-x')  ,  (jg) 

and  is  shown  in  Figure  (3a). 

In  contrast  to  the  drive  of  the  first  example,  which  was  short  and  impulsive  in  nature,  the 
drive  considered  here  is  long,  on  the  order  of  a  reverberation  time  for  the  system.  The  direct 
component  is  shown  in  Figure  (3b)  and  the  component  of  the  reverberant  response  is  shown  in 
Figure  (3c).  It  can  be  noted  that  with  the  extended  pulse,  the  system  very  quickly  settles  into  its 
modal  response.  When  all  the  response  components  are  included,  as  in  Figure  (3d),  the  system 
settles  in  to  its  modal  response  even  quicker.  The  final  figure  in  this  series.  Figure  (3e),  shows  the 
long  time  behavior  of  the  system.  Comparison  of  this  figure  with  Figure  (2e)  shows  that  at  long 
times  after  the  impulse,  the  system  response  becomes  modal  and  independent  of  the  shape  of  the 
impulse. 

The  temporal  response  of  a  system  which  is  not  stationary  in  time  is  chosen  as  the  last  of 
the  single  system  examples.  Figure  (4a)  shows  the  length  vs.  time  behavior  of  a  hypothetical 
system  which  varies  sinusoidally  from  initial  length  L  =  Lq  to  Lq  (1  ±  0.9).  The  temporal  response 
of  this  system  to  a  "short"  pulse,  equation  (17),  is  shown  in  Figure  (4b).  In  order  for  both  the 
drive  and  observation  points  to  be  within  the  system  at  its  shortest  length,  they  were  placed  at 
x'/Lq  =  0.03  and  x/Lq  =  0.07,  respectively.  Otherwise,  the  system  specifications  are  given  by 
equation  (16).  One  could  pick  any  parameter  or  combination  of  parameters  in  equation  ( 16)  to  vary 
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with  time  and  the  formalism  does  not  require  that  their  variations  be  specified  analytically;  i.e.,  they 
could  be  imperically  derived. 


MULTIPLE  SYSTEMS 

The  formalism  for  a  single  system,  presented  in  Section  3,  will  now  be  extended  to  multiple 
connected  systems.  A  set  of  systems  is  depicted  in  Figure  (5).  For  purposes  of  developing  the 
formalism,  it  is  required  that  the  individual  systems  be  spatially  one-dimensional  and  each  have  a 
propagation  velocity  which  is  not  a  function  of  the  spatial  coordinate.  Now,  however,  each  system  is 
in  an  environment  consisting  of  other  systems  with  which  it  interacts.  These  interactions  are  only 
through  the  junctions  connecting  the  interacting  systems.  In  the  spatial  formalism  [1,2],  the  set  of 
junctions  is  formulated  as  the  (square)  matrices  Ar  and  The  cross  terms  in  the  matrices,  e.g.  Arij( 
i  *  j,  describe  the  interaction  between  different  systems,  e.g.  from  the  (j)th  to  the  (i)th  at  the  r(th) 
junction;  and  the  diagonal  terms,  e.g.  Arj  j> ,  describe  the  "self-interaction."  The  former  quantities  are 
generally  termed  the  transmission  coefficients  whereas  the  latter  are  the  reflection  coefficients,  and  for 
conservative  junctions, 

^  Aaji  -  1 

a  =  r  or  q  .  (19) 

In  the  current  formalism,  the  junctions  may  also  be  characterized  with  delays  upon  reflection  and/or 
transmission.  Accordingly,  the  junction  delay  matrix  is  defined. 

Ja  =  (AajiXaji)  ,  (20) 

where  taji  is  the  delay  experienced  by  the  pulse  when  passing  from  the  i(th)  to  the  j(th)  system. 

Thus,  the  off-diagonal  elements  of  la  give  the  attenuation  and  delay  during  transmission  between 
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systems  and  the  diagonal  elements  give  the  reflection  coeffecients  and  delays.  The  delay  operator  is 
defined  to  distribute  over  the  elements  of  a  matrix,  thus,  e.g., 


A  (la)  -  (Aaji  A(xajj))  . 

Each  system,  e.g.  the  (j)th,  is  terminated  by  end  points  x,j  and  x^,  and  therefore  is  of  length 
Lj  =  |  x,j  -Xqj  |.  These  sets  of  quantities  are  formulated  in  the  vectors  xr>  xq  and  L  with 
L  =  |  xr  -  xq  |.  The  set  of  pulse  propagation  velocities  associated  with  the  systems  are  cast  into  the 
diagonal  matrices  ur  and  u4  The  drive  and  observation  points  are  also  cast  into  the  vectors  x/  and  x 
and  it  should  be  noted  that  if  a  complex  is  driven  (or  observed)  in  only  one  system,  then  the  drive 
vector  x'  (or  observation  vector  x)  has  only  one  non-zero  element  The  total  system  delays  are  cast 
into  diagonal  matrices  whose  elements,  e.g.  the  (j)th,  are  the  delay  operators  associated  with  a 
traverse  of  the  length  of  the  (j)th  system,  i.e. 


A[x(xp|xa)]=  (A  [t  (xpj  |  x a j ) ] )  8ij 


(22) 


The  reverberation  operator  corresponding  to  equation  (13b)  can  now  be  written,4 

oo 

£>a  =  n£o[A[t(xa|xp)]  A(Xp)  A  [x  (xp  |  xa)  ]  A  (Xa)  ]n  . 


(23) 


This  is  the  generalization  of  equation  (13b)  and  analogous  to  the  factor  D  in  references  1  and  2. 


With  these  definitions,  the  formalism  for  the  multiple  connected  complex  can  be  stated.  The 
operator  h  is  sought  for  which 


POL,  t)  =  h(x ,  1 1 1')  peQt',  t') 


(24) 


477ie  matrix  reverberation  operator,  equation  (23),  must  carry  an  r  or  q  subscript  since,  in  general, 
the  matrix  operators  involved  do  not  commute. 
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where  the  hjj  element  describes  the  response  in  the  (j)th  system  due  to  a  drive  in  the  (i)th  system.  The 
drive  is  again  separated  into  components  emanating  in  the  two  directions,  i.e. 


Pe  (x  >  t  )  -  L  pea(x,»  O  • 

~  a  =  r,q  ~ 

Again,  the  separation  into  direct  and  reverberant  components, 

h  (x  ,  x',  1 1 t')  =  hd  (x  ,  x',  1 1 1')  +  hrev  (x  ,  x',  1 1 1')  , 


(25) 


(26) 


with  the  direct  impulse  response  operator, 

hd  (x  ,  x',  1 1  t')ji  =  A  [tji  (xj  |xp]  =  A  [(uO'1 1  xj  -x'j  |  ]  U  (xj -x\)  5ji 
+  A  [(u^)'1 1  xj  —x'J]  U  (x'j  -  xj )  8ji  , 

(for  Xjj  >  Xqj)  .  (27) 

The  four  reverberant  impulse  response  operators  are: 

h«(&  >  t  1 1')  =  A  [ t  Oc.  I  id)  ]  A(xa)  i2>a  A  [l  (ia  |  &')  ]  , 
ha(X,i',  t  |t')  = 

A  [ T  (x  |  Xp)  ]  A(Xp)  A  [x  (xp  |  xa)  ]  A(Ta)  A  [t  (xa  |x')  ]  . 

(28) 


RESPONSE  OF  MULTIPLE  SYSTEMS:  AN  EXAMPLE 

A  simple  example  of  the  formalism  presented  in  the  foregoing  section  is  illustrated  by 
calculating  the  temporal  response  of  a  complex  consisting  of  two  connected  systems  excited  by  an 
impulse.  Two  sets  of  calculations  are  again  presented;  one  for  the  case  when  the  impulsive 
excitation  is  a  short  pulse  and  one  for  which  it  is  a  relatively  long  tone  bu^st.  The  format  of 
presentation  will  follow  the  foregoing  examples  illustrating  the  single  system. 
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The  system  parameters  used  for  the  examples  are: 


ur  =  u*t  =  u ,  (the  systems  are  isotropic), 

»  »  s: 

Xq  =  {0,  0}, 


u'1 1  =  (to,  1.4  to),  (the  temporal  "lengths"  of  the  systems), 

t|r =  T|q  =  rj  =  {0.005,  0.005},  (loss  factors), 


Ar 


(0.55  0.4 
V  0.45  0.6  ’ 


a  =  0.45  0.35 
=  0.55  0.65 


(junction  matrices), 


Tr 


(no  delays  upon  reflection 
or  transmission  between  systems), 


£  =  {0.33  Li,0},  (the  complex  is  driven  in  system  #1), 


x  =  {0.75  Li,  0},  (the  response  is  assessed  in  system  #1). 


(29) 


The  drives  for  the  two  examples  are  chosen  to  be  the  same  as  for  the  single  system 
examples,  i.e.  a  short  pulse  and  a  long  pulse  given  by  equations  (14)  and  (15)  respectively,  and 
shown  in  Figures  (2a)  and  (3a)  respectively.  Since  both  the  observation  and  the  drive  points  are 
chosen  to  be  in  the  same  system  and  are  separated  by  the  same  amount  in  both  sets  of  examples, 
the  direct  temporal  responses  for  the  two  system  examples  are  the  same  as  in  Figures  (2b)  and  (3b) 
respectively  and  these  figures  are  not  repeated.  The  reverberant  response  components  are  different 
from  the  single  system  examples.  The  p{  reverberant  component  for  the  short  pulse,  equation  (17) 
is  shown  in  Figure  (6a).  All  components  of  the  response,  the  direct  plus  the  four  reverberant  for 
the  short  pulse  excitation,  are  shown  in  Figure  (6b).  This  total  response  is  shown  again  in 
Figure  (6c)  with  the  amplitude  and  time  scales  extended  to  20  times  the  pulse  transit  time  for  the 
system  which  contains  the  drive. 
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The  same  figures  for  the  case  using  the  long  pulse  drive,  equation  (18),  are  given  in 
Figure  7;  i.e.  Figure  (7a)  shows  the  p£  reverberant  component  and  Figures  (7b)  and  (7c)  show  the 
total  response  in  the  formats  of  Figures  (6b)  and  (6c)  respectively. 

The  temporal  response  as  a  function  of  two  variables,  t  and  x,  is  shown  in  Figure  (8). 

Figure  (8a)  shows  the  spatial  and  temporal  evolution  of  the  response  in  the  driven  system  to  a  short 
pulse,  equation  (17).  The  drive  is  applied  at  t'  =  0  at  x'j/ui  i  =  0.33  x0  (in  system  number  1)  and  is 
shown  propagating  away  from  the  drive  point  and  interacting  with  the  second  system  at  the  junctions. 
The  parameters  for  the  complex  are  given  in  equation  (29)  with  the  exception  that  the  observation 
point  is  not  fixed  but  rather  used  as  a  variable.  Careful  observation  of  Figure  (8a)  reveals  response 
contributions  which  have  traversed  the  second  system  and  returned  to  the  first  system;  they  appear  as 
pulses  entering  from  the  junctions.  The  response  in  the  second,  undriven,  system  is  shown  in 
Figure  (8b). 
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Fig.  1 .  Schematic  of  a  single  one-dimensional  dynamic  system  with  terminal  coordinates  xr  and 
xq  and  reflection  coefficients  \  and  Aq  at  the  end  points.  The  drive  point  and 
observation  point  are  given  by  x'  and  x  respectively. 
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C) 

Fig.  2.  The  single  system  as  defined  by  equation  (16)  responding  to  a  "short",  unit  amplitude 

pulse  defined  by  equation  (17).  The  normalizing  time,  x0,  is  the  one-way  transit  time  for  a 
pulse  in  the  system. 

(a)  The  drive  pulse  applied  at  t  =  t '  =  0  and  defined  by  equation  (17). 

(b)  The  direct  response. 

(c)  The  p[  component  of  the  reverberant  response;  i.e.  the  response  which  results  from  the 
drive  component  emanating  from  the  point  x'  towards  the  r  junction  (c.f.  Fig.  1)  and 
which  arrives  at  the  point  x  traveling  towards  the  r  junction. 
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- 4r4~  -H  lMt - 

8.  |lo|.  I  jLRWrW.  is.  is.  ; 

(t/To)  -» 


b) 

Fig.  4.  The  single  system  with  length  varying  with  time. 

(a)  The  length  of  the  system  as  a  function  of  time. 

(b)  The  total  response,  i.e.,  direct  plus  all  four  reverberant  components  for  the 
system.  The  system  is  driven  at  (x'/Lq)  =  0.03  and  the  response  is  assessed  at 
(x/L 0)  =  0.07.  The  reflection  coefficients  are  Ar  =  Aq  =  0.95.  The  other  parameters 
are  specified  in  equation  (16). 


Pl(Xl.‘)->  (p'r)l(X,,t)-> 


Fig.  6.  The  complex  as  defined  by  equation  (29)  responding  to  a  "short”  unit  amplitude  pulse 
defined  by  equation  (17),  and  shown  in  Fig.  2a.,  applied  in  system  1.  The 
normalizing  time,  is  the  one-way  transit  time  for  a  pulse  in  system  number  1. 

(a)  The  (pr)j  component  of  the  response;  i.e.  the  response  in  system  number  1  which 
results  from  the  drive  component  emanating  towards  the  r  junction,  the  drive  being 
applied  at  the  point  xj  in  system  number  1,  and  which  arrives  at  the  point  xj 
traveling  towards  the  r  junction. 

(b)  The  total  response;  i.e.  the  direct  plus  all  four  reverberant  components. 
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Pl(Xl,  t)  -» 


Fig.  8.  (Continued) 

(b)  The  system  which  does  not  contain  the  dnve. 
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